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The exact reaction rates of beta processes for all particles at arbitrary degeneracy are derived, 
and an analytic /3-equilibrium condition Hn = i-ip + for the hot electron-positron plasma with 
nucleons is found, if the matter is transparent to neutrinos. This simple analytic formula is valid 
only if electrons are nondegenerate, which is generally satisfied in the hot electron-positron plasma. 
Therefore, it can be used to efficiently determine the steady state of the hot matter with positrons. 
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Gamma Ray Bursts (GRBs) m and core collapse su- 
pernovae(SNe) are two of the most violent events in 
our universe. Ironically, their explosion mechanisms are 
still mysterious. Recently, the central engine of GRBs is 
believed to be related to the hyperaccretion of a stellar- 
mass black hole at extremely high rates from ^ 0.01 to 



10 MqS- 



In such an accretion disk, matter is 



so dense that photons are trapped. The possible chan- 
nel for energy release is either neutrino emission which 
is mainly from the electron-positron (e^) capture on nu- 
cleons and annihilation, or outflows from the disk. 
Whatever a successful central engine is, it ejects a hot 
fireball which consists of the radiation field and baryons. 
The ratio of neutrons to protons, or equivalently, the elec- 
tron fraction, is crucial to the observed radiation from 
GRBs 0, 0, its dynamic evolution |EH and the nu- 
cleosynthesis in the disk or fireball [J, |5| . For instance, 
the inelastic collisions between neutrons and protons pro- 
duce observable multi-GeV neutrino emission 0,0- In 
addition, the two component fluid of neutrons and pro- 
tons significantly changes the fireball interaction with an 
external medium which is supposed to produce the ob- 
served electromagnetic radiation from GRBs and their 
afterglows [lo|; and the electron fraction Yg strongly af- 
fects the equation of state of the hyperaccretion disk and 
the neutrino emissions from it [ll|. 

Roughly speaking, SNe are powered by the iron core 
collapse of their progenitors. Most numerical simulations 
have shown not only the failure of the prompt shock, but 
the failure of its revival by the delayed neutrino emis- 
sion from the protoneutron star (PNS). The result (ex- 
plosion or not) sensitively depends on the input micro- 
physics, such as the electron capture, the neutrino emis- 
sion, neutrino-matter interactions, and so on (see ref.Q, 
and references therein). Without a doubt, weak inter- 
actions, especially capture and neutron decay, play a 
key role in both GRBs and SNe. During the accretion 
or collapse, these processes exhaust electrons, thus de- 
crease the degenerate pressure of electrons. Meanwhile, 
they produce neutrinos which carry the binding energy 
away. Therefore, electron capture is crucial to the for- 



mation of the bounce shock of SNe, and the resulting 
neutrino spectra strongly influence the neutrino-matter 
interactions which are energy dependent and are essential 
for collapsing simulations |l2| . 

The existence of a hot state with nucleons is the the 
common characteristics of both GRBs and SNe, as well as 
the PNS, the bounce shock and the early universe 0|. In 
these systems, the electron and positron captures are the 
two most important physical processes [lj|. The steady 
state is achieved via the following beta reactions [l5| . 
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These beta reaction rates are calculated in the previous 
studies, usually under one of three approximations: the 
nondegenerate approximation the degenerate ap- 

proximation {e.g. |l6j). and the elastic approximation in 
which there is no energy transfer to nucleons |l7l |. In 
this letter, applying the structure function formalism de- 
veloped by Reddy et al. 18] (see also J^J), I derive the 
exact reaction rates of beta processes for all particles at 
arbitrary degeneracy. In addition, I find an analytic ex- 
pression for determining the kinetic equilibrium between 
electron capture and positron capture, which is efficient 
to determine the steady state of the hot matter with 
positrons. 

Reaction rates. From Fermi's golden rule, the reac- 
tion rates of the processes read 



A = 2 



n 



d^Pi 



(2n) 



{2^f5^^\P,-Pf)\M\^T, (4) 



where Pi = {Ei,pi) denotes the four- momentum of par- 
ticle i [i = j/e/i/e, e^/e+, n,p), pi = and Pi and 
Pf are the total initial and final momentum, respec- 
tively. |M|2(= Glcos^ 6*0(1 + Sg^)) is the transition 
rate averaged over the initial spins, here Gp — 1.436 x 
10~^^ erg cm'^ is the Fermi weak interaction constant, Oq 
(sin 6*0 = 0.231) is the Cabibbo angle, and — 1-26 
is the axial- vector coupling constant. !F denotes the 
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final-states blocking factor. For instance, in reaction (1), 
•^=/e/p(l — /n), where fi is the Fermi-Dirac function of 
particle i. In this letter, we assume that the emitted 
neutrinos can escape freely from the system. Using the 
structure function formalism developed by Reddy et al. 
|l8l |. the above integrations can be simplified into only 
three dimensional ones, 
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where Si^j{qo, q) is the so-called dynamic form factor 
or structure function which characterizes the isospin re- 
sponse of the system jl^] . The expression Si^j {qo , q) is 
given by 



where 



Si^jiqo,q) = -I 7, 

1 — exp(— z) 



qo + fJ-i- fJ.j 



E' = m, 
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where /i^ and rrii are the chemical potential and the mass 
of baryons, and goi 9 denote the momentum and energy 
transfer. In Eqs. ©-Q, E^, = qo + E^, E^+ = qo + E^, 
and Ee — qo — E^, respectively. Equations ©-((TI) are 
valid for nonrelativistic and noninteracting baryons |l8l |. 
Below the nuclear density, this is a good approximation. 

Analogous to the analysis in Reddy et al. it is 

easy to obtain the previous results in the nondegener- 
ate and degenerate limits of baryons. As an illustration, 
the electron capture rate in nondegenerate limit is shown 
below, 

1 1"°° 
X^.^c^—\M\^np dE,E,p,{E,-Qff,, (12) 

where Q = m„ — nip is the mass difference between neu- 
tron and proton, iii = 2(r7iir/27r)'^/^ exp(?7i) is the num- 
ber density of neutrons and protons in the nondegenerate 
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FIG. 1: Electron capture rates on protons in SNe matter as 
a function of the baryon number density tib at different tem- 
peratures. The electron fraction Ve is taken to be 0.3. The 
solid curves show the exact results from Eq. the dashed 
curves are the elastic results from Eq. the dot-dashed 

curves are the nondegenerate results from Eq. 11211 . 



limit, and r\i = (/Zi — mi)/T is the reduced chemical po- 
tential. The above approximate rates are frequently cited 
in the literature to discuss the kinetic equilibrium for (3- 
processes and the emissivity of neutrino emissiori, even 
though its validity should be checked carefully |j, |5|, [iSj . 
In the elastic limit, rip in Eq. H12I) is replaced by 



Vpn 



= (np-nn)/(l-e(''"-'''')/^). 



(13) 



It is generally believed that electron capture rate under 
the elastic approximation in some sense introduces the 
effects of the degeneracy of baryons, so it should be more 
accurate than that under the nondegenerate approxima- 
tion. Suppose that the nuclei are dissolved completely 
into nucleons at high temperature. Figure ^ shows the 
differences between the exact electron capture rate in the 
SNe matter and the previous approximate results. The 
electron fraction Yg = {rip- — np+)/n-Q is assumed to be 
0.3, where tib = nn{^n,T)+np{^p, T) is the baryon num- 
ber density. The number density of particles at any de- 
generacy is expressed in terms of the Fermi-Dirac func- 
tions [201 . The multidimensional integrations and the 
Fermi-Dirac functions are calculated using the mixture 
of Gauss-Legendre and Gauss-Laguerre quadratures . 
From Fig. ^ it is evident that there are great differences 
between the elastic results and the exact results when 
baryons become degenerate. As shown in Eq. (|13() . the 
capture rates under the elastic approximation always de- 
crease exponentially when nucleons become degenerate, 
which is qualitatively correct in the dense nuclear mat- 
ter near /3-equilibrium. However, this conclusion is not 
correct obviously in supernova matter. The elastic ap- 
proximation is also good when /ip ~ /i„ because the fa- 
vored energy transfer is zero. The elastic approximation 
underestimates the electron capture rate, hopefully, the 
inclusion of the exact electron capture rate in the future 
collapsing simulations should be helpful for the supernova 
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explosions. 

Condition for kinetic equilibrium. In principle, if the 
matter is transparent to neutrinos which carry energy 
away, the equilibrium of the reactions can not be 

treated as a chemical equilibrium problem Suppose 
that the dynamic time scale of the system under consid- 
eration is greater than that of the reactions Q-©, the 
general condition for the kinetic equilibrium is given by 

= Ae+„(M„, Mp, Me, T) + XniHn, ^J-p, Me, T). (14) 

However, it is well known that if all the particles involved 
in the Urea processes are degenerate, such as what hap- 
pens in the interior of a cold neutron star, then the typical 
energy of emitted neutrinos is of order the temperature, 
which could be neglected compared to the Fermi energy 
of particles, i.e. , Ep + — En- So the kinetic equi- 
librium requires Xe-p — A„ which results in the so-called 
chemical equilibrium condition for the cold npe~ gas, 

Mn=Mp+Me- (15) 

Before we derive the analytic dynamical equilibrium 
condition in plasma with nucleons, we can make some 
reasonable approximations. First, the rate of neutron 
decay could be neglected before reaching the degener- 
ate limit. Second, the electrons are nondegenerate, i.e. , 
exp((i<^e — Me)/2^) > 1- Third, the energy of emitted neu- 
trinos is of order that of the captured electrons/positrons, 
i.e. , El, 2± Ee, Ef^ 2± and thus E„ ~ Ep. Under 

these approximations, X^-p = Ae+„ gives 

= Ag-p — Ag+„ OC fp{l — fn)fe+ 

/ " +1 ((g„-i;„)-(M„-M,.))/T\ 

^ /p(l-/„)/e+(e'^=/^-e(^""'''')/^). (16) 

Therefore, from Eq. (|16|l we obtain the beta-equilibrium 
condition for the plasma, 

Mn = Mp + 2Me- (17) 

It should be emphasized that during the above deriva- 
tion, it is not assumed whether the baryons are degener- 
ate or not. On the other hand, Eq. (|17|) is still valid under 
the degeneracy of baryons, which is neglected completely 
in the approximate reaction rates at the beginning. It 
is not a surprise to notice that the analytic equilibrium 
condition Eq. (|17|l can also be drawn in the completely 
nondegenerate limit. If all particles are nondegenerate, 
we have 

Ae-p OC TT-e-rip OC exp(77p -I- ?7e), (18) 
Xe+n OC ne+n„ OC exp(?7„ - ?7e). (19) 

From Eq. (fTH|l - l(TO|l . the analytic condition Eq. (fT7|l is 

also obtained hcuristically, but not very strictly. 
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FIG. 2: Upper panel: Electron fraction Fe under f3- 
equilibrium versus the baryon number density ng at differ- 
ent temperatures. The solid curves show the exact results 
from Eqs. l|Kll-||7|, the dotted curves are from the analytic 
/3-equilibrium condition Eq. 117L and the dot-dashed curves 
are from the nondegenerate approximation of baryons. Lower 
panel: The reduced chemical potentials of electrons, neutrons, 
and protons rji as a function of the baryon number density at 
different temperatures. From top to bottom, the lines corre- 
spond to the results of T = 6 x 10^ lO"', 5 x 10^°, 10"K. 



If neutrinos are trapped, the chemical equilibrium 
condition for both cold npe~ and hot npe^ gases is 
Mn + Mi'c — fJ-p + fJ-e- The chemical potential of the 
trapped neutrinos is generally assumed to be zero, thus, 
Eq. H15|l can also be understood as the chemical equi- 
librium condition for plasma with neutrino trapping 
[^. As T — > 0, the number density of trapped neutrinos 
np^ = n^^ OC T"^ — > 0. Therefore, the difference between 
Eq. (|15|l and Eq. H17|) clearly shows the effects of neutrino 
trapping. 

Validity of the analytic condition. The analytic condi- 
tion Eq. (|17|l is valid, only if electrons are nondegenerate, 
which is generally satisfied. If not, the number density 
of positrons will decrease exponentially. In the following, 
we check Eq. (|17ll by numerical calculations. As before, 
suppose that the nuclei are dissolved completely into nu- 
cleons at high temperature. Given the baryon the baryon 
number density and temperature, the electron fraction Y^. 
is determined by two equilibrium conditions. One is the 
charge neutrality n^-{^e,T) — ne+(/ie,T) = np{^p,T), 
the other is the beta equilibrium condition: the integra- 
tion equation Eq. (|14|l . or our analytic equation Eq. I|17|) . 
To check the validity of Eq. H17|) , the equilibrium Ye is cal- 
culated based on the exact reaction rates Eqs. ®-tZI), ths 
analytic condition Eq. H17|l . and the approximate rates, 
respectively. Figure |2 shows Ye and the reduced chemical 
potentials rji versus the baryon number density at differ- 
ent temperatures. It is evident that the results from the 
analytic condition are almost consistent with the exact 
results. The difference occurs only in the regime where 
electrons become degenerate. As shown in the lower 
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FIG. 3: The neutrino emissivity versus the baryon number 
density at different temperatures. The line are the same as 
those of the upper panel in Fig. |21 



neutron degeneracy is much more important than that of 
electrons, if such conditions are satisfied. 
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panel of Fig. [3 neutrons become degenerate as electrons 
do. The particles become degenerate if their reduced 
chemical potential exceeds the temperature. For elec- 
trons and neucleons, their degenerate number densities 
are estimated to be « ~ T"^ I?>'k'^ ~ 9.0 x lO-^Tfifm"^, 
and nf^ ~ (2m„,pT)i V37r^ =i 2.8 x XQ-'^Tl^ivn-^ re- 
spectively, where T\\ = T/IO^^K. In the same regime, 
there are great differences (of several orders) between our 
results and the the previous results in which the degen- 
eracy of nucleons is completely neglected. In any case. 
Fig. [3 evidently shows that the results from the analytic 
condition are much more accurate than those from the 
approximate rates in all the parameter regions. 

The advantage of having the analytic equilibrium con- 
dition at hand is obvious. For instance, we can derive 
some useful formulae under the nondegenerate approxi- 
mation. In such limit, n„/np — \ — l/Y,, = exp((^„ — 
A^p — Q)/T) = exp((2/ie — Q)/T), Therefore, we have 



2^3 

3^ 



m 



1 _^g(2Me-Q)/T- 



i20) 

The left equality is exact for the relativistic (e.g. [2l|). 
Based on Eq. (IJOJ, A'e and Ye{iie) can be determined. 
This simple result completely reproduces the previous nu- 
merical results which is corresponding to the dot-dashed 
curves in Fig. [3 At ^e/T < 1 and (2/ie -Q)/T < 1, 
Eq. H2U|I can be simplified further, 



1 (l + 0.5Q/r) 
2(l + 1.5nB/T3)' 



(21) 



A similar result to Eq. pijl was obtained previously in a 
different method 

The total neutrino emissivity under /3-equilibrium is 
shown in Fig. |3| Compared with the exact result, both 
approximate methods overestimate the rate of the neu- 
trino emission because neglecting of the degeneracy of 
particles increases the phase space for the relevant reac- 
tions. It is clearly shown in Fig. |21 that the effect of the 
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